Spintronics of a Nanoelectromechanical Shuttle 
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We consider effects of the spin degree of freedom on the nanomechanics of a single-electron tran- 
sistor (SET) containing a nanometer-sized metallic cluster suspended between two magnetic leads. 
It is shown that in such a nanoelectromechanical SET (NEM-SET) the onset of an electromechanical 
instability leading to cluster vibrations and "shuttle" transport of electrons between the leads can 
be controlled by an external magnetic field. Different stable regimes of this spintronic NEM-SET 
operation are analyzed. Two different scenarios for the onset of shuttle vibrations are found. 
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As the downsizing of electronic devices reaches the near 
molecular scale, the Coulomb forces associated with inho- 
mogeneous charge distributions produced during device 
operation become comparable with the chemical forces 
that hold the device together. Consequently, the spatial 
configuration of a device may well change dynamically 
during its operation. This inherent feature of nanoelec- 
tronics can be turned into an advantage by designing the 
devices with the mechanical degrees of freedom in mind; 
this is the scope of nanoelectromechanics and the basis 
for nanoelectromechanical systems (NEMS) [J- 

A pronounced nanoelectromechanical phenomenon — 
electron transport by means of a nanoelectromechani- 
cal "shuttle" mechanism — has recently been predicted 
to occur in certain NEMS systems as a result of a 
bias voltage-induced nanoelectromechanical instability 
0. The most suitable device for the experimental ob- 
servation of this phenomenon is a nanoelectromechanical 
single-electron transistor (NEM-SET). A NEM-SET is 
a single-electron transistor (SET) with a movable cen- 
tral island whose center-of-mass motion is confined by 
some potential. Experimental studies of NEM-SET de- 
vices can be found in Refs. They have also been 
extensively investigated theoretically [6ll7„.8„.9„ JQ|. 

Another rapidly develo ping branch of condensed mat- 
ter physics is spintronics [llj, which deals with devices 
whose functionality depends on the control and manip- 
ulation of the spin rather than the charge of electrons. 
A recent development of great interest in this context 
is the demonstration of a magnetic hybrid nanostruc- 
ture consisting of a single Cgn-molecule placed between 
two ferromagnetic electrodes [l^ ■ The possibility to ma- 
nipulate the spin of mechanically shuttled electrons in a 
molecular NEM-SET of this kind brings about an excit- 
ing opportunity to trigger nanomechanical vibrations in 
the THz range by means of a weak external magnetic 



field much smaller than the magnetic anisotropy fields 
in the leads. This Letter is devoted to exploring this 
very phenomenon. We will show that the spin-dependent 
tunneling of electrons between differently polarized leads 
results in a strong sensitivity to an external magnetic 
field of the nanoelectromechanical instability that is re- 
sponsible for the shuttling. Depending on the magnitude 
of the magnetic field two different scenarios for the onset 
of nanoelectromechanical shuttling are possible when the 
electrical field between the leads reaches its critical value. 
This demonstrates the interesting possibility to develop 
spintronics in the context of NEMS devices. 

We will consider a NEM-SET with fully spin-polarized 
magnetic leads 113|. All electrons in the left lead are as- 
sumed to have spins pointing up while in the right lead all 
spins are pointing down. The movable central island has 
a single electron energy level, which is spin-degenerate in 
the absence of a magnetic field. A symmetric bias volt- 
age is applied, resulting in an electric field £ between the 
leads. The external magnetic field B is oriented perpen- 
dicular to the direction of the magnetization in the leads. 
We consider a symmetric coupling to the leads, in which 
case the spin polarization in the leads does not contribute 
to the magnetic field on the island. 

The Hamiltonian used to describe our system. 
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has several terms. The first describes noninteracting elec- 
trons in leads {a = L,R), whose electron densities of 
states V are assumed to be energy independent. The op- 
erator a^/j(aQfc) creates (destroys) an electron with mo- 
mentum k in the lead a with the corresponding spin. 
The electrons in each lead are held at a constant elec- 
trochemical potential fiL.B. ~ TsV/2, where e < is the 
electron charge and y > is the bias voltage. Since 
the leads are fully spin-polarized the lead index a can 
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also be used as a spin index: L =] and R =[. The 
second term represents tunneling of electrons (without 
spin flip) between the island and the leads. The oper- 
ator c\{ca) creates (destroys) an electron with spin a 
in the dot. The tunneling amplitudes depend exponen- 
tially on the displacement X: Tl,r{X) = Tq exp{=FX/A}. 
The third and fourth terms describe the single electronic 
state in the dot and its couphng to the electric field £ 
and the magnetic field B. The Zeeman splitting is given 
by giisB, where g is the electronic g- factor and /is is the 
Bohr magneton. The fifth term represents the Coulomb 
repulsion between two electrons on the island and the 
last term Hose = P^/2M + describes the vi- 

brational degree of freedom associated with the center- 
of-mass motion of the island. Here M is the mass of the 
island and uq its vibration frequency. For simplicity, we 
assume that the temperature is zero. 

In the absence of an external magnetic field, electronic 
tunneling through the system described by the Hamilto- 
nian Q is blocked. Applying a bias voltage and a mag- 
netic field allows the electronic transport to be externally 
manipulated by lifting this "spin-blockade" . By deblock- 
ing electron tunneling the mechanical degree of freedom 
is also greatly influenced, the equilibrium position of the 
island becoming unstable if the rate of energy transfer 
from the electronic subsystem to the nanooscillator ex- 
ceeds a critical value. 

The coupled electronic and mechanical dynamics of the 
dot is governed by a quantum master equation for the 
corresponding reduced density operator p{t). It can be 
derived from the Liouvillc-von Neumann equation for the 
total system by projecting out the degrees of freedom 
associated with the leads and the thermal bath. The 
reduced density operator p{t) obtained in this way acts on 
the Hilbert space of the dot, which is a tensor product of 
the Hilbert space of the oscillator and the electronic space 
of the dot. The latter is spanned by the four basis vectors 
|0), I T) = c\\0), I i) = c||0), and |2) = c|c||0). In the 
electronic basis the operator p{t) can be written as a 4 x 4 
matrix whose elements are operators in vibration space. 
The diagonal elements po = (0|p|0) and p2 = (2|p|2) 
represent the density operators of the empty and doubly 
occupied oscillator correspondingly. The singly occupied 
oscillator is described by the 2x2 block pi = {p)si.s2 = 
{{si\p\s2)), where si,S2 = | T), I i)- 

In the high bias- voltage hmit (eV > U » fioj) Q 
the time evolution of the density operators po, Pi and 
P2 is determined by the coupled system of dimensionlcss 
equations of motion 

dtPo = -i[Hosc + xd,pn\-{VL{x),pi3}/2 

+ Tr,VrlJpiA/fjJ+£^po, (2) 



dtPi = -i 
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{ri,Pi}/2 



dtP2 = 



'i [Hose - xd, P2] - {Tr{x), P2} /2 
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Here all lengths are measured in units of the zero point 
oscillation amplitude xo, all energies in units of hujo and 
time in units of (^q^] x = X/xq and p = xqP/1i are 
dimensionlcss operators for the oscillator displacement 
and momentum correspondingly. Dimensionlcss electric 
(d) and magnetic {h) fields are defined by 



d — e£ / [Mlj^xq) , h = gpsB / {ThjJq) 



(5) 



The tunneling of electrons is described by the di- 
mensionlcss parameters Ta{x) = 2TrVT^{x + d)/{'huJo), 
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damping of vibrations is introduced via the simplest form 
of the damping Liouvillian, 



-ij [x, {p, .}] /2 - 7(n^„ + 1/2) [x, [x, .]] /2 , (6) 



where 7 <^ 1 is a dimensionlcss dissipation rate and 



1] is the Bose distribution function. 



One can formally derive Eq. JHJ by weakly connecting 
the oscillator to an Ohmic heat bath by adding the terms 

Hbath = Y,q^^qblbq and Hosc-bath = XY,q9q[b\ + bq) 

to the Hamiltonian and then eliminate the bath vari- 
ables in the Born-Markov approximation . 

The mechanical degree of freedom alone is described by 
the density operator = po + T^spi + P2, which is con- 
veniently analyzed in the Wigner representation defined 
by 
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d^e-^^Hx + ^/2\p\x~^/2) (7) 



for a density operator p. 

In an experimentally relevant regime where the 
electromechanical coupling parameter 



r/ = d/\ - 1/A < 1 



(8) 



is small, one can find the stationary solution of the sys- 
tem of Eqs. (|2I3I4|I perturbatively in terms of the small 
parameters 77 and 7. 

Rescaling the phase variables X = x/A, P = p/X 
and changing to the polar coordinates X — Asimp, 
P = Acosip, one obtains the steady-state equation for 
the Wigner function W+ 



Wp_^_ of the oscillator as 



{-T+iX)C/2 + C^}W+ 
{Tjdp+T^iX)C/2}Wo-2 



(9) 

r+(X)CWx-T/2. 



Equation ^ couples to the steady state equation for the 
vector-function W = [W0-2, W^i-T' W^-c W+c, W0+2] , 



where Wo±2 



Pll 
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= Wp^^ ± Wp^^. Here r±(X) = TniX) ± r^CX), 
C^=-f [dpP + (1/2A2)9|.] and C = cosh [i\-'^dp] - 1. 

One can show that in the leading order approximation 
the Wigner function W+ {A, ip) is (/^-independent and de- 
termined by 



OaA [f{A) + D{A) Oa] W+{A) = , 



(10) 



where f{A) = (A/2) [7 - 7?/3o(A)], A(3o{A) = 
— {dip / n) cos 6*0*^2 (V') ^ and D{A) > is of sec- 
ond order in rj and 7. The function Gq*^2(¥') is deter- 
mined by the system of differential equations: 



(11) 



D^Gfl^ - r_(X)Gr2 - - T+iX) , (12) 

i?^GL°] = 2hGfl^, D^ = 2d^ + r+{X) (13) 

It follows from Eq. pU|) that the Wigner function W+ 
has the form W+{A) « Z^^ exp |- /g^ S^} and is 
peaked for amplitudes Am determined by the conditions 
/{Am) = and /'{Am) > 0. In the vicinity of these am- 
plitudes W+ {A) can be approximated by a narrow Gaus- 
sian of variance = D{Am)/ f (Am)- The behavior 
of the stationary solution W+(yl) is determined by the 
structure of the positive definite function f3o{A), which 
is bounded, has only one maximum and decreases mono- 
tonically for large A. One can show that if /i < the 
function /3o(A) bas its only maximum at A = 0, while if 
h > he it has a minimum there. In the vicinity of A = 0, 
one finds that Pq{A) = ''jthr/v + C(^^)i where 
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From this equation it follows that the electromechani- 
cal properties for a low-transparency junction (with a 
resistance in the Gohm range) is sensitive to very weak 
magnetic fields of order 1-10 Oe. Such weak fields have 
a negligible effect on the internal magnetization of the 
leads, which is why this effect was not considered here. 

The presence of various stationary regimes can be il- 
lustrated by a "phase diagram" in the {d, /i)-plane. Fig- 
ure ^ shows three domains that correspond to three dif- 
ferent types of behavior of the nanomechanical oscillator. 
In the "vibronic" domain (v), defined by the condition 
77/7 < 1/ [max/3o(A)] (h), the NEM-SET system is stable 
with respect to mechanical displacements of the island 
from its equilibrium position. The "shuttle" domain (s) 
is determined by the condition 7 < jthr and corresponds 
to developed island vibrations of a classical nature. The 
third domain is the "mixed" domain (m). It appears be- 
cause the V- and s-regimes become unstable at different 
combinations of electric (d) and mag netic {h) fields if h 
exceeds a critical value he (cf . 0, Q| ) . 

Hence, while the shuttle regime is unstable below the 
line V in Fig. ^ the vibronic regime becomes unstable 




FIG. 1: "Phase diagram" (for r+(0) = 0.1) in the electric field 
(d)- magnetic field (h) plane (see Eq. ((SJ) showing domains of 
different nanoelectromechanical behavior of the studied NEM- 
SET device. In the vibronic domain (v) only the vibrational 
ground state of the device is stable, in the shuttle domain (s) 
the stable state corresponds to developed island vibrations 
and in the mixed domain (m) both states are locally stable 
with probabilities that are equal along the dashed line p. 



above the line s. Between these lines (in the m-domain) 
both states can be stable. The oscillator "bounces" be- 
tween the V— and s-regimes due to random electric forces 
caused by stochastic variations of the grain charge asso- 
ciated with tunneling events. 

The transition time between the two locally sta- 
ble regimes of the m-domain is given by t^^s = 
ljq^ exp {Sy^s/v)^ where Sy^sid^h) ^ 1. Since Sy^s 
Ss^y and 77 <^ 1, the difference between the switching 
rates Ty—,s and t„^s can be exponentially large. This 
implies that the probabilities for the system to be in one 
or the other of the two regimes can be very different. 
The line p in Fig. ^ corresponds to Sy^s — Sy^s and 
therefore to equal rates for the transitions w — > s and 
s — > D. Below this line, the probability for the system 
to be in the w-regime is exponentially larger than for it 
to be in the s-regime, while above the dominance of the 
s-regime is exponentially large. Due to the smallness of 
the electromechanical coupling, 77 <C 1, the transition be- 
tween the two regimes is very sharp. Hence the change 
of vibration regime can be regarded as a "phase transi- 
tion" . Such a transition will manifest itself if the varia- 
tion of the external fields is adiabatic on the time scale 
of max{T<;<^^}. One can expect enhanced low- frequency 
noise, u < t^^^, around the line p as a hallmark of the 
transition. 

In the opposite non- adiabatic limit, cither the s- or the 
w-regime is "frozen in" in the mixed domain after cross- 
ing the line p. Thus if one starts in the u-domain the 
u-regime persists until the system crosses the line s, and 
if one starts from the s-domain the s-regime persists until 
the system crosses the line v. Hence, one observes a hys- 
teretic behavior of the non-adiabatic shuttle transition. 

As one can see from Fig. ^ there are two different see- 
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narios for the onset of shuttle vibrations. If one crosses 
over from the v- to the s-domain when h < i.e. avoid- 
ing the mixed domain, the onset is soft. In this sce- 
nario the vibration amphtude grows gradually from zero 
to some finite value after crossing the border line (Fig.|2Il. 
If /i > he, on the other hand, the onset is hard. In this 
case the vibration amplitude shows a step at the transi- 
tion point (Fig. ISJ, which corresponds to crossing either 
the p- or the s-line depending on whether the transition 
is adiabatic or not. 
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FIG. 2: Steady-state amplitude Ac of shuttle vibrations and 
current I through the NEM-SET device for parameters cor- 
responding to a "soft" transition between the vibronic and 
shuttle domains of Fig.[Tl(r+(0) = 0.1 and d/{^X) = 14.3). 
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FIG. 3: Hysteretic behavior of the steady-state shuttling am- 
plitude Ac and current / through the NEM-SET device for 
parameters corresponding to a "hard" transition between the 
vibronic and shuttle domains of Fig. (r+(0) = 0.1 and 
d/(7A) = 11.1). 

One can show that the expression for the steady state 
current through the system is 

1 = ^ j j dXdPTL{X)[W+ + Wfi^2 + Wi-^]/2. (15) 

Typical plots of the current in the cases of soft and hard 
transitions are shown in Figs. [3 and 13 respectively. 

In conclusion, we have considered "shuttle" phenom- 
ena in a nanoelectromechanical single-electron transistor 



(NEM-SET) system consisting of a metallic island sus- 
pended between spin-polarized leads. We have shown, 
that a coupling between the transport of spin-polarized 
electrons and the center-of-mass motion of the island al- 
lows us to control the dynamics of the mechanical de- 
gree of freedom of the island by an external magnetic 
field. Different stable operating regimes of the magnetic 
NEM-SET were found and transitions between them in- 
duced by varying the electric and magnetic fields were 
analyzed. We have hence demonstrated that magnetic- 
field-controlled spin effects can lead to a very rich behav- 
ior of nanomechanical systems. 

Although we have considered fully spin-polarized leads 
and assumed a symmetric set-up, the overall qualitative 
picture does not change if these conditions are relaxed. 
A partial polarization, allowing a shuttle instability even 
at zero external magnetic field, makes the system sensi- 
tive to the magnetic field in a more narrow interval of 
electric fields. An asymmetric coupling to the leads in- 
duces an additional magnetic field along the direction of 
their spin-polarization. This field reduces the rate of en- 
ergy pumping into the oscillator, moving the domains in 
Fig. n to higher values of electric and magnetic fields. 
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